Orbital ordering in undoped manganites via a generalized Peierls instability 
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We study the ground state orbital ordering of LaMnOz, at weak electron-phonon coupling, when 
the spin state is A-type antiferromagnet. We determine the orbital ordering by extending to our 
Jahn- Teller system a recently developed Peierls instability framework for the Holstein model^. By 
using two-dimensional dynamic response functions corresponding to a mixed Jahn- Teller mode, we 
establish that the Q2 mode determines the orbital order. 
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I. INTRODUCTION 

Undoped manganites like LaMnO^ are the parent 
systems for the colossal magnetoresistive materials. It 
is well known that orbital ordering occurs around 
780 K resulting in a C-type orbital structure with 
two kinds of orbitals alternating on adjacent sites in 
the xy plane while like orbitals are stacked in the z 
direction 2 -. As the temperature is further lowered to 
140 K, an A-type spin antiferromagnetic order sets in 
wherein the spins are ferromagnetically aligned in the 
xy plane with the spin coupling in the z direction be- 
ing antiferromagnetio 3 ". To explain the observed or- 
der several studies have been reported. These stud- 
ies fall into two broad classes based on the dominant 
cause for the observed order. One class corresponds 
to electron-electron (Coulombic) interaction^^ 7 , being 
the main cause while the other class treats the coop- 
erative Jahn- Teller (JT) interactio n 8 ' 9 ' 10 ! 11 as the more 
important one. Lin and Millis 12 have made a quanti- 
tative analysis of the effects of both interactions, gen- 
erally concluding that both pieces of physics are impor- 
tant, but with many subtleties. There is further contro- 
versy about the strength of the electron-phonon inter- 
action with extended X-ray absorption fine structure 13 
and pulsed neutron diffraction 14 measurements point- 
ing to strong interaction while some electron microscopy 
measurement a 15 ' 16 have inferred weak coupling in the 
charge ordered phases. In that regime optical measure- 
ments often infer small electronic gapsi 7 - and measure- 
ments of nonlinear transport 1 ^ have been interpreted as 
due to sliding motion of a density wave. 

Without addressing ah initio the issues of the quan- 
titative strength of the interactions it is worth under- 
standing how in principle a weak coupling theory might 
possibly work. The notion of JT is a molecular one, and 
the linear splitting of levels by a local distortion a use- 
ful principle only if the induced gap is much larger than 
the bandwidth (which it is not). Nonetheless, oxides 
are generally viewed as a template for strong interac- 
tion physics, both of the electronic and phononic variety. 
In this paper, we step back from the complexities of the 
full many-body theories to point out that the canoni- 



cal model for LaMnO^ has a weak- coupling generalized 
Peierls instability that reproduces qualitatively the or- 
dering observed. One advantage of the simplification in- 
troduced by our approach is that we can study effects of 
adiabaticity that turn out to enter logarithmically in the 
ratio of electronic bandwidth to phonon frequency. 

Our observation follows straightforwardly from assum- 
ing A-type antiferromagnetic ordering. On account of 
strong Hund's coupling, the transport is restricted to 
spin polarized electrons in two dimensions only, where 
furthermore the bands are strongly nested. The proxim- 
ity to a nesting instability allows us to employ the weak- 
coupling framework developed earlieri 1 ^ and analyze the 
orbital ordering by using a generalized Peierls instability 
approach. However, as compared to the one-dimensional 
Peierls charge density wave (CDW) approach, our higher 
dimensional orbital density wave (ODW) analysis is more 
complicated on account of there being two e g orbitals 
(with inter-orbital hopping) and two response functions 
corresponding to the JT Q2 and Q3 distortions. The 
consequences of a nesting instability on the orbital or- 
dering in LaMnOz were first discussed by Yarlagadda 
and Mitra3°- and later qualitatively by Efremov and 
Khomskii2i. 

In this paper, we study the Peierls instability condition 
by extending the recently developed reliable condition in- 
volving the dynamic susceptibility^ to a mixed JT mode. 
We find that Q2 Jahn- Teller distortion, as observed ex- 
perimentally, preempts other JT normal mode distortions 
at all values of adiabaticity and temperature. Further- 
more, the condition of instability (i.e., functional depen- 
dence of critical coupling on adiabaticity) is qualitatively 
similar to that of the one-dimensional single-orbital Hol- 
stein model. Lastly, we also find that mean-field approx- 
imation (in spite of being crude) and static Peierls insta- 
bility condition (albeit erroneous) indicate that Q2 mode 
rather than Q3 mode determines the orbital order. 



II. MODEL HAMILTONIAN 

We will now consider manganite systems with two e g 
orbitals per site and ignore spin. The Hamiltonian con- 
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while the phononic part of the Hamiltonian is given by 



FIG. 1: Fermi seas corresponding to the eigenenergies A^ 



sists of the kinetic term, the ionic term, and the electron- 
ion interaction term. The kinetic term in momentum 
space is given by 



Bp*, 



(1) 



where Bt = (b\p, b^p) with b\ and 62 corresponding 
to the destruction operators for electrons with the or- 
thonormal wavefunctions ipx 2 -y 2 an d ip3z 2 -r 2 respec- 
tively. Furthermore, T is a hermitian matrix with T^i = 
— 1.5t[cosp x +cosp y ], T 2 ,2 = —0M[cosp x +cosp y ], and 
Ti,2 = 0.bV3t{cosp x — cospy]. The eigenvalues of the 
kinetic energy are given by A^ = —t [cos p x + cosp y + 
(— l) n ^cos 2 p x + cos 2 p y — cosp x cospy] with n — 1,2. 
The Fermi sea corresponding to the lower eigenenergy 
value A2 is given by the union of the region — it/ 2 < 
k x < 7r/2 (with all values of k y allowed) and the region 
— 7r/2 < k y < ir/2 (with all values of k x allowed) as 
shown by the shaded region (both dark and light) in Fig. 
[TJ Whereas the Fermi sea corresponding to the higher 
eigenenergy value A J;' is given by the intersection of the re- 
gion — 7r/2 < k x < ir/2 and the region —tt/2 < k y < ir/2, 
i.e., only the dark shaded region in Fig. [T] Since the 
number of electrons is equal to the number of sites, the 
total area occupied by both Fermi seas is equal to the 
area of the Brillouin zone (47r 2 ). Furthermore, the Fermi 

surface corresponds to = 0. 

The electron-phonon interaction term is given by 
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+Q^[ j bx j -b\ j b 2j )], (2) 



Ho 



J '=2,3 



(3) 



where + fj- = \j2MuJaQij. Since we are interested in 
understanding orbital order, our Hamiltonian does not 
contain breathing mode distortions. 



III. PEIERLS INSTABILITY 

In this section, to understand orbital ordering at weak 
electron-phonon coupling, we consider the Peierls insta- 
bility condition by using the dynamic susceptibility in- 
stead of the static one (see Appendix A for a justifica- 
tion). The cooperative Jahn- Teller effect requires com- 
patible distortions on adjacent sites which implies that 
the ordering wavevector in two-dimensions is given by 
Q = (tt, 7t). We expand the free energy to quadratic or- 
der in the relevant degree of freedom (i.e., density n of 
electrons in an appropriate occupied orbital) as follows: 



F = 



q=±Q 



+ gu}on-f[(fl_-) + (f H )\ 
(4) 



where fyj + /f • = \/2MujqQ $j . Here Q<j> is the dominant 
mode defined as = Q 3 cos(20) + Q 2 sin(2<-/>) where 
only orbitals tp x 2_ y 2 cos((j))+ip z 2 sin(</>) or their orthonor- 
mal orbital states —ip x 2_ y 2 sin((/>) + tp z 2 cos(<^) are occu- 
pied. The order parameter corresponding to phonons is 
given by (/^) = |(/ g)|e je . Using reflection symmetry, 

we first note that x<p{Q,u) = X<f>(—Q> w )j n Q = n -Qi 
(Uq) = (W, and ( f y = (fl_ ). For g > 0, free 
energy minimum occurs at = n. Minimizing F, with 
respect to |(/^)|, yields |(/^)| = gn^. Thus, we get 



F = - 



1 + 2g 2 LU ReX4>(Q> u} ) 



Rex<p(Q,u) 
On defining the effective susceptibility as 



eff 



(5) 



(6) 



1 + 2g 2 uj Rex4> ' 
the Peierls instability condition is given by 

1 + 2g 2 u Bex4,(Q,uo) = 0, (7) 

and leads to the divergence of x% (Qi w o)- We take u> = 

m xl^iQ'U) because u is the natural frequency for 
lattice distortion. A better explanation for choosing u> = 
luq is given in Appendix A. 

We need to determine at what value of <fi one gets the 
largest value of Hex^iQ , wo) • Then one can determine 
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which normal mode gives the lowest value of g = g c sat- 
isfying the Peierls instability condition. Note that, as 
the rotational angle <j> (for < < 7r/2) is varied, all 
possible normal modes are spanned starting from Q3 at 
(\> = to Q2 at <fi = 7r/4 and then to —Q3 at <fi = ir/2. 
The dynamic susceptibility is given by 



X<f> 



(<7»=E 



KnlpQMm 2 \(0\ P QM\n)\ 2 



,(8) 



where 

P QM = £ [b\ u ^ k - b\ % J^ cos(20) 

k 

+ \bK Xt + bK J> lE | sin(2<7>). (9) 

L lk+q 2/c 2k+q IfcJ V ^' W 

Then, after some algebra, one gets 



Rex^^o) = J! 

k,a,/3 



X fe+f k+q' 



LU 



\k \k+q 



COS 



'k + q+ 6 k + ( a + P> 



2(j> 



,(10) 



where a = 1,2; p = 1,2; (ei^cf) = (6^,6^) ■ M, 
M is the diagonalizing matrix for the kinetic matrix 
T with M M = sm(0g/2), M 2 , 2 = -sin(^/2), and 
Mi. 2 = cos(0g/2). It is interesting to note that, for 
symmetric wavevectors q — (q, q) , there is no coupling 
between the density operators corresponding to Q 2 and 
Q3 modes because the inter-orbital hopping Ti : 2 = 
0.5^/3t[cosp x — cospy] is asymmetric with respect to in- 
terchange of momenta p x and p y . Thus for q = (q, q), we 
obtain 

xAq^o) = X3(<7>o)cos 2 (20) +X2(g,wo)sin 2 (20), (11) 
where X2,3 correspond to JT modes Q2,3 
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FIG. 2: Plot of X2,3{Q,Q) as a function of scaled temperature 
aT at at = 1.0 eV and Q = Q = (7T,7r). 



X2,3(Q, 0) vary logarithmically with ksT/t for t/ksT > 2 
and thus have the form 



Re[-*X2,3(Q,0)] = m 2 , 3 \n{t/k B T) + k 2 - 



(12) 



We find that m 2 « 12.6 (m 3 w 4.2) and k 2 w 18.3 (k 3 w 
18.5) with the ratio 77x2/7713 taking the expected value 
3. Thus it appears that Q2 mode is likely to dictate the 
orbital ordering. 



B. Dynamic Instability Case 

While both the static Peierls instability condition and 
the mean-field energy analysis (see Appendix C) depend 
only on the polaron size parameter (g 2 LOo/t), here for the 
dynamical Peierls instability condition [of Eq. (J7J)] there 
are two relevant parameters - namely adiabaticity pa- 
rameter t/uiQ and electron-phonon coupling g. We find 
that for any value of the adiabaticity parameter t/ujo the 
maximum value of Rex<f>{Q, c^o) occurs at cf> = 7r/4 which 
corresponds to Q2 mode. In Fig. [3l using Eq. (jTUJ) , a 



A. Static Instability Case 

Now, although the static Peierls instability condition 
1 + 2g 2 uiQX<t>(Q, 0) = erroneously predicts instability 
even for vanishing values of g, it can still help identify 
which normal mode produces the Jahn- Teller instability. 
We will first present results for the static susceptibilities 
X2, 3 (Q, 0). From the plot of X2,a(Q> °) (shown in Fig. |SJ) 
as a function of scaled temperature aT (with a being a 
scaling parameter and hopping term at set equal to 1.0 
eV) we see that they diverge logarithmically as T — > 
with %2 diverging faster than \3- At K, both X2(Q, 0) 
and X3(Q,0) produce a divergence because of the fact 
that Xi +< ^ = — A§ and that the Fermi energy is zero. 
Furthermore, the ratio X2(Q, 0)/Xs(Q> 0) = 3 at K (see 
Appendix B for details). As can be seen from Fig. [21 
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FIG. 3: Plot of Rex</> as a function of the adiabaticity t/ujo 
for various values of 4>. 4> = 0° (45°) corresponds to X3 ixs)' 
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FIG. 4: Plot of the critical coupling g c as a function of adia- 
baticity i/a>o for the susceptibilities \2 an d X3- 

variation of Rex<p(Q, ujq) (at K) is plotted for a few rep- 
resentative values of cj) — 0,7r/12,7r/6,7r/4. The curves 
for Rex w /i2(Q, Wo) and Rex 7r / 6 (Q, lu ) (in Fig. [3]) ver- 
ify Eq. pip . Furthermore, we also found numerically 
that Rex<p(Q, ojq) [given by Eq. (fPU)) ] is symmetric about 
<f> = 7t/4 - a fact that follows from Eq. (fTT|) . 

Quite strikingly, all the Rex<p(Q, u>o) vary logarithmi- 
cally with the adiabaticity t/uio for t/u>o > 2 and have 
the form 

^ e [-tX(k(Qi w o)] = m<p ln(t/o;o) + k^. (13) 

We find that m^/A w 12.6 (mo ~ 4.2) and K,y 4 w 25.5 
(/to ~ 20.9). Interestingly, the slopes in Eq. (flU|) 
are the same as those in Eq. (|T2l . The ratio of the 
slopes m n u/mo = 3 as expected from the fact that 

X2(Q, 0)/X3[Qi 0) — 3 at K. Furthermore, this logarith- 
mic dependence is quite like that for the Holstein model. 
Using the dynamic Peierls instability condition, similar 
to the Holstein model case, we are lead to an instability 
condition of the form lu = a\te~ a2t l 9 w ° where are 
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FIG. 5: Plot of the susceptibilities Rex2,3(Q, Wo) as a function 
of the scaled temperature aT for values of scaled hopping 
at = 1.0 eV and adiabaticity t/coo = 5, 100. 



constants. We also calculated the critical value of the 
electron-phonon coupling g c at which the instability oc- 
curs if only Q2 mode or only Q3 mode is excited. We find 
that the value of g c increases monotonically with the adi- 
abaticity parameter (similar to the Holstein model) and 
that, as expected, the g c value is the smallest for Q2 dis- 
tortion (as can be seen in Fig. [3J at any value of t/ojQ. 

We have also studied the temperature dependence of 
the dynamical susceptibilities (as shown in Fig. [5]) and 
find that at low temperatures the curves are constant 
with the extant of the constant region increasing as t/ujQ 
decreases. Such a behavior is consistent with the ex- 
pectation that RexfiiQ , t^o) is constant over the region 
ksT << luq. Furthermore, at higher temperatures the 
susceptibilities for various adiabaticities merge. For in- 
stance, when aT attains a value of around 300 K, curves 
for t/u>o = 100 and 00 merge (as can be seen from Figs. [2] 
andO; and for aT around 2000 K, curves for t/uo = 100 
and 5 merge. The high temperature behavior too is un- 
derstandable because one expects the effect of non-zero 
value of loq to vanish when k B T >> ujq. At the Jahn- 
Teller orbital ordering temperature of 780 K and for real- 
istic values of both t and ujq (i.e., for 0.15 eV < t < 0.38 
eV-2 2 . and for 0.06 eV < iv < 0.07 eV 23 ), range of the 
critical coupling [as obtained from Eq. ([7]) and Fig. [5] is 
0.2 < g c < 0.28. For instance, at T = 780 K, t = 0.2 eV 
[and hence a = 5 in Fig. (J5J)], and uj = 0.07 eV, we get 
g c w 0.21. Lastly, for k B T » uj but k B T/t < 0.5, the 
curves display a logarithmic dependence on ksT/t which 
is in tune with the logarithmic dependence on ujq ft of the 
susceptibility of the Holstein model when uio/t < 0.5 (see 
Ref. B)- 



IV. CONCLUSIONS 

We will now discuss the general features of the orbital- 
ordering instability and compare it with the Peierls in- 
stability in the Holstein model. For the Holstein model, 
at K, the mean-field approximation gives a gap A of 
the formal 

A = 8te' 7rt/g2uJ0 . (14) 
In our case as well, we find that the gap is given by 

A 2) 3«4 3 ie-<3t/s 2 "o ! ( 15 ) 

where d 2 ' 3 are constants and A2,3 ~ <7 2 wo|c2,3| with 
C2.3 being amplitudes of orbital density waves defined 
in Appendix C. It should however be noted that, when 
A/ujq << 1, mean-field gives erroneous results. For in- 
stance, it predicts a gap even when the electron-phonon 
coupling g is small. Although, mean-field approximation 
is inaccurate at the transition, it can still help us figure 
out which of the two JT modes is dominant. As shown 
in Appendix C, mean- field correctly shows that Q2 mode 
prevails over Q3 mode. 
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Next, in the Holstein modeli*^, at K and t/ujQ > 2, 
the actual instability condition is given by 

lo = 8te-^ B ' U0 . (16) 

For our JT system too, at UbT « uo {k B T >> lo ) and 
when t/u>o > 2 (t/ksT > 2), the instability is of the form 

Wo (7^r) = ai te-»^», (17) 

and thus, like the Holstein model, has an essential sin- 
gularity at g = 0. We also note that one cannot get the 
correct Peierls instability condition by the approximation 
P|/(2M) + KQ\/2 w KQ 2 /2 for the normal mode dis- 
tortion even when t/u>o is large. This is because, when 
P?/(2M) = 0, the double commutator for the distortion 
Qtf, becomes zero, 

Q^ = -[[Q^H],H} = 0, (18) 

which implies that phase transition always occurs! 

In summary, we observe that owing to the one- 
dimensional like Fermi surface at zero doping in man- 
ganites (as shown in Fig. 1), there are strong similari- 
ties of the above mentioned nature between our JT sys- 
tem and the one-dimensional Holstein model. The one- 
dimensionality of our manganite system is a result of the 
flatness of the Fermi surface [as can be seen, for instance, 
from Eq. jB7}]. When t/[max{u; , k B T, A }] > 2, we 

find that the susceptibility Re[— tx4 > {Q 1 <^o)] varies log- 
arithmically with respect to t/[max{uJo, k B T, A^}] and 
has the general form 

Re[-t%^,((5, Wo)] = rn^ \n(t/[max{uj , ^k B T, A^}]) + k^, 

with 7 « 1.77 and both and being given by Eq. 
(fT5|) . Using this logarithmic relation and the generalized 
Peierls instability condition of Eq. ([7]), one obtains the 
explicit form of the instability condition. 

In conclusion, we have studied orbital ordering for 
the ground state of the undoped manganite systems in 
the weak electron-phonon coupling regime gu>o/t < 1. 
We employ the generalized dynamic Peierls instabil- 
ity condition 1 + 2g 2 ujQY{,ex<t>{Q, loq) = to figure out 
which normal mode or combination of normal modes 
causes the instability. It is also important to note that 
the dynamic Peierls instability condition does not suf- 
fer from the problem of predicting CDW instability at 
vanishingly small electron-phonon coupling (i.e., g — ► 
0) as does the usual static Peierls instability condition 
[1 + 2g 2 uj Q x4Q,0) = 0]. We find that Q 2 Jahn-Teller 
distortion produces the first instability and thus pre- 
empts other normal mode distortions. Thus the two- 
dimensional orbital ordering, in the ferromagnetic planes 
of the observed A-type antiferromagnetic state, is gov- 
erned by the Q2 JT mode being cooperatively excited in 
the system. Hence, we find that the experimentally ob- 
served order can be explained even without considering 
electron-electron interactions. 



Before we close, a few general discussions are in order. 
Above the magnetic transition temperature Tjv, where 
orbital structure does not change much, transport is per- 
mitted in the third direction and the Fermi surface for 
three dimensions should be considered. Then, although 
the bands are not flat, we still have the nesting condi- 
tion X-y = — A§ for Q — (7r,7r,7r) and hence the static 
susceptibilities will diverge. However, the experimental 
ordering wavevector is (tt, 7r,0) and not (7r,7r,7r). To get 
the observed ordering one will have to incorporate addi- 
tional physics such as octahedral tilting. Next, at non- 
zero temperatures below Tjv, hopping in the third direc- 
tion is small but non-zero owing to the non-saturation 
in A-type antiferromagnetic order. Then flatness (one- 
dimensionality) of the Fermi surface would be lost. How- 
ever, hopping in the third direction increases with tem- 
perature and the situation is different from that men- 
tioned in Ref. |2^ where, since the hopping in the trans- 
verse direction decreases with increasing temperature, re- 
entrant behavior could occur. Lastly, electron-electron 
interactions can have an effect on the nesting conditions 
as pointed out by Kugel, Sboychakov, and Khomskii 27 . 
These authors find that electron-electron interactions 
lead to the occurrence of nesting at a density of less than 
an electron per site. However, in this work, Luttinger's 
theorem is violated [see Fig. 5(e) in Ref. H3| and implica- 
tions of that should be investigated for non-Fermi liquid 
behavior. If, indeed in a full-fledged calculation, beyond 
the Hubbard I approximation, nesting (with flat Fermi 
surface) occurs at a lower density, then a corresponding 
ODW instability condition should be re-analyzed for such 
a situation. 
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APPENDIX A 

We shall give a heuristic justification for the use of 
dynamic susceptibility in the Peierls instability condition 

l + 2g 2 Lu Rex^Q,Lu ) = 0. (Al) 

Let Q^q be the dominant normal mode distortion oper- 
ator (at wavevector Q) in the Fourier transformed space. 
We know that the double time derivative of the operator 

Q&q is § iven h y 

Q^ = -[[Q^H},H}. (A2) 
Then on taking matrix elements we get 

(*i|Q*3l*o) = -(E^ - E^fi^lQ,^), (A3) 
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where $„ is an eigenstate with n phonons all of which 
are in the state Q. When aA.. = (E$ 1 — E< in ) 2 < 0, in- 
stability occurs for transition from |$o) to |$i) provided 
that ($1 IQ^qI^q) 7^ 0. Now, at weak electron-phonon 
couplings (i.e., when guj^/t < 1) 

E<s> 1 — E$ = ojci + ReS0(<5, tt>o) 

= cj + 5 2 ^Rex0(Q, w ), (A4) 

where is the self-energy corresponding to mode Q^q- 
Thus, when 

cu 2 eff = wg[l + 2.g 2 w o Rex0(Q,c;o)] = 0, (A5) 

CDW instability occurs. The above instability condi- 
tion is exact up to second-order in perturbation theory. 
A more detailed and rigorous derivation of the dynamic 
Peierls instability condition is given in Ref. [l|. 



APPENDIX B 

We will show analytically that X2(Q, 0)/x3(Q, 0) = 3 
at K. Understanding the susceptibilities is complicated 
because the eigenstates [corresponding to the eigenval- 
ues Xi 2 ] are a linear combination of the states ipk,x 2 -y 2 
and ipk,3z 2 -r 2 with coefficients that are a function of the 
wavevector k. More precisely, the eigenvectors for X^ 2 
are given by (c 1 - ,c^) = (b^?, ■ M, where M is the di- 
agonalizing matrix for the kinetic matrix T with M^i = 
sin(0 £ /2), M 2i2 = -sin(%/2), and M u = cos(%/2). 
Now, from the the kinetic matrix T, we get 



cos(6>p-) = 



0.5[cospa; + cospy] 



■J cos 2 p x + cos 2 p y — cospz cos p y 



and 



sin 



O.Sv^cosp.j — cos p y ] 

J COS 2 p x + COS 2 p y — COS Pa; COS p y 



(Bl) 



(B2) 



In the expressions for X2,3(Q,0) given below 

x 2 (Q,o) 



E 

k,a,0 



(S -) 

X k+Q k+Q' 



k+Q 



sin 



A* - A 



h+Q + °k + ( a + P)* 



, (B3) 



and 



x 3 (Q,o) = 



E 

k,a,/3 



x k k' 



(S a 

v k+Q k+Q 1 



k+Q 



xl - X 



h+Q + h + ( a + P)* 



, (B4) 



because A 2 +( ^ = — X\ and since on the Fermi surface (FS) 
A j = 0, the following term diverges 



(clVA - (c 2f -c 2 A 

x k k' v k+Q k+Q 1 



x\-x 



k+C 



(B5) 



Furthermore, because X* = — A 2 and since A 2 = on 
the FS, the following term also diverges 



( C 2t c 2 ) - ( c y -c 1 - -) 

x k k' v k+Q k+Q 1 



A 2 



V fc+Q 



(B6) 



Then 



X2(Q,o)/x 3 (Q,o) = 



cos 2 


~ k + Q+ e k~ 




2 


FS 


sin 2 


~ e k + Q+ 6 k 




2 


FS 


sin 2 ( 




3, 


cos 2 ( 


e k)FS 



(B7) 



where use has been made of the fact that the FS is 
flat and one-dimensional like and that on the FS either 
k x = ±7r/2 or k y = ±7r/2. 



APPENDIX C: MEAN-FIELD CDW ANALYSIS 

Assuming that the total wavefunction of the system is 
separable into a phononic part and an electronic part, 
after averaging the Hamiltonian over the phononic coor- 
dinates, we get the following effective Hamiltonian (with 
details given in Ref. Il9l ): 

^=E b p T % 

-2g 2 u J2 [(bljhj + OiJ !>\M, + &2>y) 
i 

+ (blihj-bljhMibii-btjhj) 
+g 2 cu Y / (b\ J b2 1 + b^hjf + (bl^ - b^) 2 , (CI) 



where < .. > implies averaging over the relevant coordi- 
nates which here are electronic. 

Based on the arguments that wavevector Q determines 
the orbital ordering in two-dimensions (as discussed in 
Sec. Ill), we compute the ground state energy using 
mean-field when only either Q 2 mode or Q% mode gets 
excited cooperatively in the system. The order param- 
eters are given by (b\jb 2 j + b\^hj) — c 2 cos(Q ■ Rj) and 
(bljhj - b\jb 2 j) = c 3 cos(Q • Rj) with -1 < c 2j3 < 1 
and Rj being the position vector. Here it should be 
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FIG. 6: (a) Dependence of dimensionless ground state energy 
per site (E/t) on dimensionless polaronic energy (g 2 uio/t) for 
cooperative Q2 and Q3 modes; (b) variation of coefficients 
£2,3 of ODW order parameters for Q2 and Q3 distortions as 
a function of q 2 ujo /t. 

pointed out that the order parameter (b\j 623 + ^jy^ij) 
corresponds to the density difference of electrons in the 
two orbitals ipx = (i J x 2 -y 2 ~ V'32 2 -r 2 )/v / 2 and tpy = 
— (^ x 2_ v 2 + 7/j 3z 2_ r 2)/v / 2 (as described in Ref. IToh . 



The unit cell needed to compute the ground state en- 
ergy consists of two adjacent sites with the Brillouin zone 
being given by — it < (k x + k y ) < it and — it < 
(k x — k y ) < 7T. We diagonalize a 4 x 4 matrix at each 
momentum and integrate the lowest two eigenenergies 
over the Brillouin zone to obtain the ground state en- 
ergy. The results of our calculations are shown in Fig. 
[5] From Fig. HJa) we see that the ground state energy 
corresponds to the Q2 mode with the difference in energy 
between the Q 2 only state and the Q 3 only state peaking 
at intermediate values of the dimensionless polaronic en- 
ergy {g 2 ujQ 1 1) . For zero values and infinite values of the 
polaronic energy both modes yield the same energy be- 
cause zero value implies no phononic coupling effect while 
infinite value corresponds to localized polarons. Thus for 
large values of the polaronic energy, the ground state en- 
ergy is only slightly smaller than the polaronic energy. 
Furthermore, from Fig. [6](b) we also see that, as the po- 
laronic energy increases, the values of increase and 
become unity around g 2 u>o/t ~ 2 implying that for the 
Q3 (Q2) mode ip x 2_ y 2 (ipx) orbital is occupied fully at 
one site with the tp 3z 2_ r 2 (ipy) orbital being fully occu- 
pied at the adjacent sites. 



1 S. Datta and S. Yarlagadda, Phys. Rev. B 75, 035124 
(2007). 

2 Y. Murakami, J. P. Hill, D. Gibbs, M. Blume, I. Koyama, 
M. Tanaka, H. Kawata, T. Arima, Y. Tokura, K. Hirota, 
and Y. Endoh, Phys. Rev. Lett. 81, 582 (1998). 

3 E. O. Wollan and W. C. Koehler, Phys. Rev. 100, 545 
(1955). 

4 K. I. Kugel and D. I. Khomskii, Sov. Phys. JETP 37, 725 
(1973). 

5 S. Ishihara, J. Inoue, and S. Maekawa, Phys. Rev. B 55, 
8280 (1997). 

6 L. Sheng and D. N. Sheng, Int. J. Mod. Phys. B 13, 1397 
(1999). 

7 S. Okamoto, S. Ishihara, and S. Maekawa, Phys. Rev. B 
65, 144403 (2002). 

8 A. J. Millis, Phys. Rev. B 53, 8434 (1996). 

9 T. Hotta, S. Yunoki, M. Mayr, and E. Dagotto, Phys. Rev. 
B 60, R15009 (1999). 

10 P. B. Allen and V. Perebeinos, Phys. Rev. B 60, 10747 

(1999) . 

11 Z. Popovic and S. Satpathy, Phys. Rev. Lett. 84, 1603 

(2000) . 

12 C.Lin and A.J Millis, Phys. Rev. B 78 174419 (2008) 

13 A. Lanzara, N. L. Saini, M. Brunelli, F. Natali, A. Bian- 
coni, P. G. Radaelli, and S.-W. Cheong, Phys. Rev. Lett. 
81, 878 (1998). 

14 Despina Louca, T. Egami, E. L. Brosha, H. Roder, and A. 
R. Bishop, Phys. Rev. B 56, R8475 (1997). 

15 J. C. Loudon, S. Cox, A. J. Williams, J. P. Attfield, P. 
B. Littlewood, P. A. Midgley, and N. D. Mathur, Phys. 



Rev.Lett. 94, 097202 (2005). 

16 S. Cox, E. Rosten, J. C. Chapman, S. Kos, M. J. Caldern, 
D.-J. Kang, P. B. Littlewood, P. A. Midgley, and N. D. 
Mathur, Phys. Rev. B 73, 132401 (2006). 

17 K. H. Kim, S. Lee, T.W. Noh, and S.-W. Cheong, Phys. 
Rev. Lett. 88 167204, (2002). 

18 S. Cox, J. Singleton, R.D. McDonald, A. Migliori, and PB. 
Littlewood, Nature Materials 7, 25 (2008). 

19 S. Yarlagadda, Int. J. Mod. Phys . B 15, 3 529 (2001). 

S. Yarlagadda and M. Mitra, arXiv:cond-ma t/0310350| 
(unpublished) . 

21 D. V. Efremov and D. I. Khomskii, Phys. Rev. B 72, 
012402 (2005). 

22 A. E. Bocquet, T. Mizokawa, T. Saitoh, H. Namatame, and 
A. Fujimori, Phys. Rev. B46, 3771 (1992); T. Arima, Y. 
Tokura, and J. B. Torrance, Phys. Rev. B48, 17006 (1993); 
T. Saitoh, A. E. Bocquet, T. Mizokawa, H. Namatame, A. 
Fujimori, M. Abbate, Y. Takeda, and M. Takano, Phys. 
Rev. B51, 13942 (1995). 

23 M. N. Iliev, M. V. Abrashev, H.-G. Lee, V. N. Popov, Y. 
Y. Sun, C. Thomsen, R. L. Meng, and C. W. Chu, Phys. 
Rev. B57, 2872 (1998). 

24 R. H. McKenzie, C. J. Hamer, and D. W. Murray, Phys. 
Rev. B 53, 9676 (1996). 

25 R. J. Bursill, R. H. McKenzie, and C. J. Hamer, Phys. 
Rev. Lett. 80, 5607 (1998). 

26 J. H. Wei, D. Hou, X. R. Wang. larXiv:0811.2019l 

27 K. I. Kugel, A. O. Sboychakov, D. I. Khomskii, J. Super- 
cond. Nov. Magn. 22, 147 (2009). 



